Back

Psychological Medicine

Cambridge University Press (CUP)

Preprints posted in the last 7 days, ranked by how well they match Psychological Medicine's content profile, based on 74 papers previously published here. The average preprint has a 0.07% match score for this journal, so anything above that is already an above-average fit.

1
Interoceptive accuracy and attention across multimorbidity classes: A latent class analysis

Mulder, J.; Boeker, C. M.; Smit, A. K.; Kiefte-de Jong, J. C.

2026-06-09 public and global health 10.64898/2026.06.08.26355147 medRxiv
Top 0.1%
12.5%
Show abstract

Background Multimorbidity is increasingly prevalent, and associated with worse clinical and psychosocial burdens. Interoception, the brain's ability to sense and interpret internal bodily signals, may contribute to multimorbidity, through its link with health behaviors, stress regulation, and mental health. This study examines whether self-reported interoceptive accuracy and attention is associated with multimorbidity, by identifying multimorbid subgroups and their interoceptive profiles. Methods Morbidity classes were identified through latent class analyses in two Dutch survey datasets, focusing on depression and alexithymia (DA-dataset; N = 671) and lifestyle factors (L-dataset; N = 1022). Linear regression analyses were used to assess interoceptive accuracy and attention (by the Interoceptive Accuracy Scale and Interoceptive Attention Scale respectively) among different subgroups. Results Multimorbid subgroups were characterized by older age, low socioeconomic position, and elevated physical, psychological, and behavioral problems. Multimorbid classes exhibited lower interoceptive accuracy (DA-dataset: B = -1.14, 95% CI = [-2.89, 0.62]; L-dataset: B = -2.36, 95% CI = [-3.83, -0.89]) and higher attention (DA-dataset: B = 3.62, 95% CI = [0.97, 6.27]; L-dataset: B = 1.07, 95% CI = [-1.42, 3.56]) compared to healthier classes. Conclusion Multimorbid populations demonstrated lower interoceptive accuracy and higher interoceptive attention. This highlights the psychosocial complexity of multimorbid populations which may impact their self-management and health behavior. These findings underscore the need to expand treatments to include psychosocial domains for multimorbid patients.

2
Mortality in people with attention-deficit/hyperactivity disorder (ADHD): Examining how risk is embodied in a pooling of two prospective cohort studies

Li, H.; Ford, T.; Warrier, V.; Bell, S.; Batty, G. D.

2026-06-09 epidemiology 10.64898/2026.06.08.26355148 medRxiv
Top 0.1%
12.4%
Show abstract

Background. Nascent findings suggest that people with attention-deficit/hyperactivity disorder (ADHD) experience higher rates of mortality. To date, study samples have been insufficiently well-characterized to examine the mechanisms via which this neurodevelopmental condition elevates mortality risk. Methods. We used data from the 2007 and 2011 waves of the US National Health Interview Survey, a general population-based cohort study comprising 52097 adults (28675 women) aged 18 years or older at baseline. ADHD diagnosis and an array of demographic, socioeconomic, lifestyle, and co-morbidity (somatic and psychiatric) covariates were self-reported. Findings. At baseline, compared with unaffected individuals, participants with ADHD were more likely to be socioeconomically disadvantaged, smoke cigarettes, consume alcohol, and report symptoms of psychological distress. A median 7.75 years of mortality surveillance (range: 7.25-12.25) gave rise to 6597 deaths from all-causes. After adjustment for age, sex, ethnicity, and survey year, ADHD was associated with a markedly elevated risk of death (hazard ratio [95% confidence interval]: 1.58 [1.20-2.09]). Statistical adjustment for socioeconomic circumstances (11% attenuation), physical co-morbidities (15%), and lifestyle factors (17%) had only a modest impact on the ADHD-death gradient, with the greatest explanatory power apparent for symptoms of depression and anxiety (58%). The magnitude of the association of ADHD with mortality was commensurate to that for several well-established risk factors such as poverty (1.66 [1.55-1.78]), hypertension (1.41 [1.32-1.51]), and diabetes (1.71 [1.59-1.85]) but somewhat lower than cigarette smoking (2.51 [2.29-2.76]) after controlling for age, sex, ethnicity, and survey year. Associations between ADHD and cause-specific mortality from cardiovascular disease, cancer, and chronic respiratory disease were inconclusive. Interpretation. In the present study, the influence of ADHD on total mortality appears to be largely embodied via a series of malleable characteristics, particularly mental illness. If confirmed elsewhere, these results raise the possibility that risk factor modification via standard pharmacological and behavioral interventions could help reduce rates of premature mortality in this patient group. Funding. This paper received no direct funding. GDB is supported by the UK Medical Research Council (MR/P023444/1) and the US National Institute on Aging (1R56AG052519-01, 1R01AG052519-01A1).

3
Exploring the role of binge eating in the association between ADHD and BMI: A twin study

YOU, Y.; McAdams, T.; Oginni, O.; Liu, C.; Herle, M.; Zavos, H.

2026-06-05 psychiatry and clinical psychology 10.64898/2026.05.28.26354354 medRxiv
Top 0.1%
9.7%
Show abstract

Objective: ADHD has been associated with obesity indicators, including BMI, across the lifespan. A possible mechanism linking ADHD and BMI is binge eating. Previous research has found associations between ADHD, binge eating and BMI. However, the role of genetic and environmental influences on these associations remains unclear. Method: We utilized data from the Twins Early Development Study (TEDS), comprising 3,675 monozygotic and 7,063 dizygotic twin pairs. ADHD symptoms in childhood and adolescence were assessed using parent-reported questionnaires. Adult ADHD symptoms were measured using both self-report and parent-report questionnaires. Phenotypic mediation models examined whether binge eating mediated the association between ADHD and BMI, without controlling for genetic confounding. Subsequently, the etiological architecture underlying the associations among the three traits across childhood, adolescence, and adulthood were investigated by incorporating genetic and environmental influences into the models. Results: Binge eating significantly mediated the association between ADHD symptoms and BMI in both adolescence and adulthood. However, these mediation effects were no longer present once genetic and environmental influences were incorporated into the models. The best-fitting model in childhood, adolescence and adulthood was Cholesky decomposition models, where covariance between traits was explained by shared aetiology. Conclusions: This twin study reveals shared liability across ADHD, binge eating, and BMI. The mediating role of binge eating in the relationship between ADHD symptoms and BMI was largely confounded by shared genetic influences. Intervention strategies could focus more on common underlying behavioural and self-regulatory mechanisms across these traits, as well as placing more emphasis on symptom patterns within families.

4
Developmental Associations Linking Childhood Trauma and Early Cannabis Use to Adolescent DNA Methylation and Psychotic-Like Experiences

Trotta, G.; Liu, Z.; Austin-Zimmerman, I.; Spinazzola, E.; Sideli, L.; Aas, M.; Rodriguez, V.; Li, Z.; Leung, B. M.; Li, Q.; Zhang, S.; Sham, P. C.; Vassos, E.; Bentall, R.; Walker, E. M.; Dempster, E.; Murray, R.; Di Forti, M.; Alameda, L.; Wong, C. C. Y.

2026-06-10 psychiatry and clinical psychology 10.64898/2026.06.09.26355257 medRxiv
Top 0.1%
9.6%
Show abstract

Background. Psychotic-like experiences (PLEs) index early risk for psychotic disorders and are consistently associated with childhood trauma, yet underlying biological mechanisms remain poorly understood. DNA methylation (DNAm) may capture the biological embedding of early adversity, while adolescent exposures such as cannabis use may modify these processes. We examined epigenome-wide associations of childhood trauma and PLEs, tested the moderating role of early cannabis use, and evaluated DNAm as a potential mediator. Methods. We analysed data from the Avon Longitudinal Study of Parents and Children (ALSPAC), a UK population-based birth cohort. Childhood trauma was assessed prospectively and retrospectively. Epigenome-wide DNAm was measured in peripheral blood at ~17 years using the Illumina 450K array, and PLEs were assessed at 18 using a structured interview. Epigenome-wide association studies were conducted for trauma-DNAm and DNAm-PLEs associations in the final sample (n = 1,457), adjusting for demographic, biological, and technical covariates. Differentially methylated regions (DMRs) were identified using DMRff, followed by functional enrichment analyses. Cannabis use at 15.5 was modelled as a moderator with multiple imputation for missing data. Mediation was tested using the Divide-Aggregate Composite-null Test (DACT). Results. Childhood trauma was associated with widespread DNAm differences, primarily at the regional level, with enrichment in pathways related to cellular stress responses. In contrast, DNAm associated with PLEs was more limited and implicated loci involved in epigenetic regulatory processes. These signatures were largely distinct, and there was no evidence supporting mediation after multiple testing correction. Incorporating cannabis use altered the pattern and extent of DNAm associations, with stronger and more significant signals observed at both CpG and regional levels, although these did not translate into evidence of mediation. Conclusion. Childhood trauma and PLEs show distinct DNAm signatures in adolescence, with trauma-related DNAm reflecting broad stress-related processes and PLE-associated DNAm implicating regulatory mechanisms. We found little evidence that DNAm mediates the trauma-PLE association. Instead, adolescent exposures, particularly cannabis use, may distinctly influence trauma-related epigenetic variation with limited detectable downstream effects on PLEs. These findings support a context-dependent model of epigenetic risk and highlight the need for larger longitudinal studies to clarify causal pathways linking early adversity to psychosis.

5
Neuroimaging Summary Scores Predict Trajectories of Psychotic-Like Experiences in Youth

Cooper, R. E.; Sahasrabudhe, R.; Glahn, D. C.; Jalbrzikowski, M.

2026-06-04 psychiatry and clinical psychology 10.64898/2026.06.03.26354754 medRxiv
Top 0.2%
8.2%
Show abstract

Objective. Persistent, distressing psychotic-like experiences (PLEs) are associated with neurobiological alterations and increased psychosis risk. We combined individual-level neuroimaging measures with effect sizes from large neuroimaging studies to create a summary score ('Psychosis Neuroscore') reflecting neuroanatomic liability for psychosis, and examined its ability to predict PLE trajectories in young adolescents. Method. Using latent growth mixture models, we estimated PLE trajectories from four annual visits of the Adolescent Brain Cognitive Development Study (N=9584, ages 9-10 at baseline). Using baseline T1-weighted and diffusion-weighted imaging data, we calculated Psychosis Neuroscores, as well as Neuroscores for two psychiatric disorders with late adolescent/adult onset (Major Depressive Disorder, Bipolar Disorder). We compared Psychosis Neuroscores to i) other psychiatric Neuroscores, ii) modifiable risk factors, and iii) established risk factors in predicting trajectory membership. Results. We identified four trajectories of distressing PLEs: Persistent Elevated (N=1,968, 21%), Gradual Decreasing (N=3,424, 36%), Rapid Decreasing (N=1,593, 17%) and Low/No Distress (N=2,599, 27%). Adolescents with Persistent Elevated PLEs had significantly higher Multimodal (combined T1 and diffusion-weighted) and T1-weighted Psychosis Neuroscores than all other trajectories (Odds Ratios [ORs] 1.27-1.34,pFDR<.01). Bipolar Disorder Neuroscores showed a similar pattern (ORs 1.16-1.23,pFDR<.01). Psychosis Neuroscores showed comparable associations with established risk factors in predicting trajectory membership, but smaller associations than modifiable risk factors, including screen time, physical activity, and sleep disturbances. Conclusion. Psychosis Neuroscores differentiate youth with persistent PLEs from those with decreasing, remitting or low PLEs, demonstrating their potential utility for early risk stratification. Integration with established risk factors may enhance psychosis risk prediction in youth.

6
Tune In or Take the Stage? A Randomized Controlled Trial Comparing After-School Music and Theatre Training with Neuroimaging Outcomes for Youth

Jamey, K.; Herschel, E.; Noel, C.; Villanueva, J.; Reyes, M.; Hsu, E.; Ilari, B.; Mack, W.; Luo, S.; Habibi, A.

2026-06-05 public and global health 10.64898/2026.06.03.26354844 medRxiv
Top 0.4%
4.4%
Show abstract

Introduction: While growing evidence suggests that music training supports child development, few long-term randomized controlled trials (RCTs) have rigorously tested these claims. Moreover, it remains unclear whether the benefits are confined to music-specific domains or extend to higher-order cognitive functions such as inhibitory control (IC), a core executive function associated with long-term outcomes in academic achievement, career success, socio-emotional health, and physical well-being. This paper presents the protocol for the Extracurricular Activity and Child Early Learning and Development (EXCEL) trial, an RCT designed to assess the feasibility of a long-term music training program focusing on the brain and behavioral correlates of IC. Methods: A total of 126 children, aged 6 to 8 years and residing in neighborhoods with limited resources in Los Angeles, were individually randomized to either a music (intervention) or theatre (active control) after-school program. Both programs were delivered over 24 months by established community arts organizations. Eligibility criteria included: average intellectual functioning, no major medical or psychiatric conditions, and MRI eligibility. Children with prior formal music training exceeding six months or severe hearing impairment were excluded. Before the intervention began, all participants completed baseline behavioral and neuroimaging assessments. The primary trial aim was to assess the effects of extended music training, relative to theatre training, on changes in measures of IC (i.e., Go/No-Go task and delayed gratification) and related neural functional activation. A secondary interim aim of the trial was to evaluate the feasibility of conducting a long-term RCT of music education in a first cohort, measured by participant retention, adherence to the program, willingness to continue at the 12-month mark, and fidelity. Progress: Recruitment, screening, baseline testing, randomization, and program enrollment began in August 2022, and after-school programming began in October 2022. The randomized interventions and all data for the first cohort (N = 42) have been collected. Intervention and active control programs for a second cohort are ongoing and will end in Fall 2026. Discussion: This paper reports the EXCEL trial protocol and provides feasibility estimates for implementing a long-term randomized controlled trial of music training in real-world, community-based settings with children. While similar neuroimaging RCTs are currently underway in Europe, the EXCEL trial is among the first in the United States to integrate longitudinal neuroimaging with arts intervention. Findings will inform the viability of scaling such programs and contribute to our understanding of how sustained music engagement may influence the development of inhibitory control circuitry in childhood.

7
Trans-ancestry genome-wide association meta-analysis of antidepressant response to selective serotonin reuptake inhibitors in clinical studies of depression

Hu, K.; Lo, C. W. H.; Awasthi, S.; Pain, O.; Singh, M.; Ahn, Y.; Aitchison, K. J.; Baune, B. T.; Biernacka, J. M.; Bondolfi, G.; Carrillo-Roa, T.; Choi, H.; Czamara, D.; Domschke, K.; Fabbri, C.; Hamilton, S. P.; Ising, M.; Jang, Y.; Kato, M.; Kim, D. K.; Kim, D.; Lee, B.-C.; Lewis, G.; Lim, S.-W.; Liu, Y.-L.; Myung, W.; Perroud, N.; Serretti, A.; Tsai, S.-J.; Uher, R.; Weinshilboum, R.; Won, H.-H.; Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium, ; Ripke, S.; Coleman, J.; Lewis, C. M.

2026-06-04 genetic and genomic medicine 10.64898/2026.06.03.26354703 medRxiv
Top 0.6%
3.5%
Show abstract

Antidepressants are widely prescribed for major depressive disorder, yet only one-third of patients achieve remission after initial treatment. Previous genome-wide association studies (GWAS) of clinically assessed antidepressant response combined multiple antidepressant classes, potentially obscuring class-specific effects. This study focused on selective serotonin reuptake inhibitors (SSRIs), often first-line due to better tolerability. Data from 15 cohorts across four ancestries were integrated: European (N = 3887; 11 studies), East Asian (N = 1068; 4), African (N = 277; 1), and Admixed American (N = 250; 1). GWAS of non-remission and percentage improvement were conducted within cohorts, followed by ancestry-specific meta-analyses and trans-ancestry meta-regression. Single nucleotide polymorphism (SNP)-based heritability was estimated in European samples. Polygenic scores were used for leave-one-out prediction and to assess shared genetic architecture with psychiatric traits. Gene-level and gene-set enrichment analyses were also performed. No genome-wide significant variants were identified for either outcome in any ancestry-specific or trans-ancestry analyses. However, trans-ancestry meta-regression yielded eight independent loci with suggestive associations (p < 1 x 10-5) for non-remission and 17 for percentage improvement. Gene-set analyses revealed nominal enrichment of the serotonergic synapse pathway for non-remission. SNP-based heritability estimates were not significantly different from zero for either outcome. Better SSRI response was nominally associated with lower genetic predisposition to major depressive disorder, post-traumatic stress disorder, and schizophrenia. This study represents the largest trans-ancestry GWAS of SSRI response, highlighting emerging biological signals. Limited power emphasises the need for larger and ancestrally diverse cohorts to better characterise the genetic architecture of antidepressant response.

8
Shared epigenetic regulation acting on neuroimmune pathways contributes to the comorbidity between generalized anxiety disorder and COVID-19

Karaca, S.; Cabrera Mendoza, B.; He, J.; Qiu, D.; Davtian, D.; Lacobelle, A.; Nunez, Y. Z.; Krystal, J. H.; Pietrzak, R. H.; Gelernter, J.; Polimanti, R.

2026-06-04 genetic and genomic medicine 10.64898/2026.06.03.26354830 medRxiv
Top 0.7%
2.7%
Show abstract

Background: The biological mechanisms linking generalized anxiety disorder (GAD) and COVID-19 remain poorly understood, despite substantial evidence of their comorbidity. To address this gap, we examined genetic and epigenetic factors underlying their co-occurrence. Methods: In a multi-ancestry sample of 893 participants, we conducted genome-wide and epigenome-wide analyses of GAD and COVID-19 severity. Integrating large-scale genome-wide datasets and information regarding methylation quantitative trait loci, complementary analytic approaches were used to identify regional methylation patterns, assess genetically regulated DNA methylation in blood and brain tissue, and evaluate causal loci shared between GAD and COVID-19. Results: GAD was associated with epigenome-wide significant variation in loci involved in chromatin regulation and synaptic signaling. Conversely, COVID-19-related epigenetic signals were enriched in immune-inflammatory and host-response pathways. Mild COVID-19 was epigenetically related to endothelial-inflammatory signals, while severe COVID-19 was linked to epigenetic changes implicated in myeloid and thrombo-inflammatory pathways. Epigenetic signals shared between GAD and COVID-19 implicated processes related to stress adaptation and tissue homeostasis. Genetically informed analyses identified 60 shared loci, including MAPT, ZFP57, and FBXL18, indicating pleiotropy between GAD and COVID-19 in genetically regulated DNA methylation variation. Brain-specific analyses further highlighted convergence in additional loci (i.e., MICB and HLA-DPB1), suggesting neuroimmune mechanisms underlying GAD-COVID-19 shared methylation patterns. Conclusions: These findings support that GAD and COVID-19 share epigenetic and genetic architecture involving pathways related to vascular integrity, immune function, and cellular adaptation, highlighting a potential neuroimmune basis for their co-occurrence.

9
Diagnosis provision by young people's mental health services: a comparison with epidemiological data

Lewis, S. J.; Meehan, A. J.; Akiba, M.; Arseneault, L.; Byford, S.; Caspi, A.; Clark, B. R.; Downs, J.; Ford, T. J.; Fisher, H. L.; Koenen, K. C.; Moffitt, T. E.; Newbury, J. B.; Odgers, C. L.; Pritchard, M.; Simonoff, E.; Danese, A.

2026-06-05 psychiatry and clinical psychology 10.64898/2026.05.28.26354156 medRxiv
Top 0.8%
2.0%
Show abstract

Background Little is known about the provision of diagnoses to young people with mental health disorders. We investigated diagnosis provision by NHS mental health services, focusing on 17-year-olds in South London between 2009-2024, and compared with estimated disorder prevalence. Methods To examine diagnosis provision in the population, we extracted diagnosis data from records of the NHS mental healthcare provider serving South London, using the Maudsley Biomedical Research Centre Clinical Record Interactive Search application; we then compared these data with the corresponding population size, obtained from the Office for National Statistics. To assess diagnosis provision in those with mental health disorders, we compared diagnosis data with the number of young people estimated to have met criteria for a disorder, derived from epidemiological interview data collected in the Environmental Risk (E-Risk) Longitudinal Twin Study and weighted according to characteristics of 17-year-old South Londoners. To assess diagnosis provision in those with mental health disorders within health services, we compared diagnosis data with the number estimated to have met criteria for a disorder and used any health service for their mental health, again derived from weighted E-Risk Study data. Findings Of 17-year-olds from South London in 2009-2024, 4.0% (n=8,958/223,404) had a diagnosis in mental health records during the previous year. This diagnosis provision covered <1 in 16 of those estimated to have had a mental health disorder, and <1 in 4 of those estimated to have also used health services. Diagnosis provision was lower in girls than boys and in young people with Black/Asian/Mixed/Other ethnicity than those with White ethnicity, in those estimated to have had a mental health disorder and used health services. Interpretation These findings demonstrate gaps and biases in mental health diagnosis provision for young people, including within health services, and reveal the imperative need to strengthen young people's mental healthcare.

10
Characterizing artificial intelligence (AI) psychosis in a large academic medical setting: evidence of the new clinical phenomenon and the vulnerability of those in early phases of psychosis

Bergson, Z.; Vassall, S. G.; Wright, A.; McCoy, A. B.; Schafer, K. M.; Achee, M. C.; Sheffield, J. M.

2026-06-08 public and global health 10.64898/2026.06.04.26354939 medRxiv
Top 0.9%
1.8%
Show abstract

Background: Concerns about "AI psychosis" have swirled in the media since ChatGPT's release, but few systematic analyses exist. We therefore conducted an electronic health record (EHR) analysis to identify the frequency, clinical characteristics, and quality of AI interactions in patients experiencing psychosis treated in a medical center. Methods: AI keywords (e.g., ChatGPT, AI) were used to search Vanderbilt University Medical Center's EHR from 12/1/2022-4/1/2026. Records were discarded if they were not AI-related or if the primary diagnosis did not include psychosis. Three raters read notes to determine if a patient was experiencing AI psychosis and classified the interactions using 4 a-priori categories (Catalyst, Amplifier, Co-Author, Object) formulated to explain how AI-related negative outcomes emerge. Findings: 73 patients met our criteria. 28 patients were rated as experiencing AI psychosis, 17 had neutral interactions, and 28 expressed delusional content related to AI without documented evidence of conversational AI use. ChatGPT was the matching keyword for 53.6% patients experiencing AI psychosis. The majority of AI psychosis cases were documented after ChatGPT's "4o" model was released in May 2024. Notably, the AI Psychosis group had significantly more patients experiencing a first psychotic episode (60.7%) compared to the other two groups. Amplifier was the most common (64.3%) qualitative rating in the AI Psychosis group. Interpretation: "AI psychosis" is an infrequent but real phenomenon observed in clinical practice. Most affected patients were experiencing their first psychotic episode and presented with AI psychosis following the release of the more sycophantic GPT-4o. Among the affected patients, AI most often exacerbated an existing condition by reinforcing distorted ideas.

11
Neuroanatomical dimensions in recent-onset depression: clinical profiles, inflammatory markers, and proteomic ageing

Lalousis, P. A.; Moles, L.; Antoniades, M.; Xiao, W.; Couch, A. C. M.; Erus, G.; Thokachichu, P.; Srinivasan, D.; Fan, Y.; Woodham, R. D.; Arnone, D.; Arnott, S. R.; Chen, T.; Choi, K. S.; Fatt, C. C.; Frey, B. N.; Frokjaer, V. G.; Ganz, M.; Godlewska, B. R.; Hassel, S.; Ho, K.; McIntosh, A. M.; Qin, K.; Rotzinger, S.; Sacchet, M. D.; Savitz, J.; Shou, H.; Stolicyn, A.; Strigo, I.; Strother, S. C.; Tosun, D.; Victor, T. A.; Wei, D.; Wise, T.; Zahn, R.; Anderson, I. M.; Deakin, J. F. W.; Craighead, W. E.; Dunlop, B. W.; Elliott, R.; Gong, Q.; Gotlib, I. H.; Harmer, C. J.; Kennedy, S. H.; Knudse

2026-06-04 psychiatry and clinical psychology 10.64898/2026.06.01.26354320 medRxiv
Top 1.0%
1.7%
Show abstract

Background: Major depressive disorder (MDD) is clinically heterogeneous, hindering identification of reproducible biomarkers. Using a semi-supervised machine learning approach, HYDRA, we previously identified two neuroanatomical dimensions from structural MRI in medication-free MDD from COORDINATE-MDD consortium. These dimensions (D1, D2) showed differential responses to selective serotonin reuptake inhibitor (SSRI) antidepressants and placebo. External replication in UK Biobank linked D2, characterized by widespread subtle neuroanatomical reductions, to an immuno-metabolic profile. Here, we examined whether these dimensions are detectable early in the course of illness. Methods: We applied the pre-trained model to structural MRI data from the multisite PRONIA cohort, comprising individuals with recent-onset depression (ROD; n = 377; mean age 25.8 years, SD 6.0; 51.3% female) and healthy controls (n = 267; mean age 25.5 years, SD 6.4; 61.0% female). Participants were assigned to clusters (C1, C2) corresponding to the previously identified dimensions (D1, D2). Clusters were compared on clinical symptom profiles, peripheral inflammatory markers, and in a subset (n = 107), proteomic ageing indices. Results: Two neuroanatomical clusters were identified in PRONIA. C1 (n = 265) showed higher negative symptom severity and elevated interleukin-2 levels. C2 (n = 140) was associated with higher residual proteomic age. Overall depressive symptom severity did not differ significantly between clusters. Conclusions: Neuroanatomical dimensions of MDD are reproducible and detectable at illness onset. Associations with negative symptom severity, inflammatory signalling, and proteomic ageing suggest these dimensions capture biologically meaningful heterogeneity early in depression. These findings support a biologically informed framework for stratified treatment approaches in MDD.

12
Longitudinal brain structural changes during clozapine treatment: associations with neuroreceptor architecture and clinical response

King, B.; Cannon, D.; Crossley, N. A.; Valderrama, A. G.; Hallahan, B.; Jung, W. H.; Kempton, M. J.; Kim, S.; Lawrence, A. J.; MacCabe, J. H.; McDonald, C.; Mena, C.; Nakajima, S.; Papale, A.; Raminfard, S.; Sarpal, D.; Sim, H.; Tronchin, G.; Tuominen, L.; Kim, E.; Egerton, A.

2026-06-10 psychiatry and clinical psychology 10.64898/2026.06.06.26354980 medRxiv
Top 1%
1.3%
Show abstract

In treatment-resistant schizophrenia, clozapine treatment has been associated with longitudinal reductions in subcortical volumes, ventricular enlargement, and widespread cortical thinning. However, it is unknown how these structural changes relate to clozapines pharmacological profile and clinical efficacy. We combined five longitudinal datasets with MRI acquired before and on average 5 months after clozapine initiation in 143 individuals to quantify brain structural changes and their association with normative maps relating to neuroreceptor architecture and physiological systems, and improvement in symptom severity. Clozapine treatment was associated with grey matter volume reductions across multiple subcortical regions (including the amygdala, hippocampus, thalamus, caudate, putamen and nucleus accumbens), increases in pallidal volume, ventricular enlargement, and widespread cortical thinning. Cortical regions showing the greatest magnitude of thinning corresponded to areas with higher normative densities of serotonergic 5-HT1A, 5-HT2A and 5-HT4 receptors. Changes in subcortical volume or cortical thickness during clozapine treatment were not associated with changes in total or positive symptom severity. In addition, baseline subcortical volume, cortical thickness, or gyrification prior to starting clozapine did not predict subsequent symptom improvement. Cortical thinning may partly reflect clozapines activity at serotonergic receptors, which have been implicated in cortical network stabilisation and neuroplasticity, however structural remodelling during clozapine treatment may reflect a process independent from its clinical efficacy in improving core symptoms of psychosis.

13
Direct and mediated effects (DME) SLCMA: a novel method for life course modelling with time-varying covariates

Beer, S.; Simpkin, A. J.; Eldeeb, S. Y.; Zar, H. J.; Stein, D. J.; Dunn, E. C.; Smith, A. D. A. C.

2026-06-06 epidemiology 10.64898/2026.05.29.26354427 medRxiv
Top 1%
1.2%
Show abstract

Background: In prospective cohort studies, where an exposure is collected repeatedly, interest often lies in determining whether the timing of that exposure has a differential effect on a later outcome. The Structured Life Course Modeling Approach (SLCMA), where users select between temporal hypotheses of exposure specified a priori, provides one way to analyse such longitudinal data. However, few studies using SLCMA consider the effect of time-varying covariates (TVC) which may impact associations. Methods: We present a modified version of the SLCMA - called direct and mediated effects (DME)-SLCMA - which corrects for TVC. We first develop the DME-SLCMA method, test it through simulation, and apply it to psychosocial data from the Drakenstein Child Health Study (DCHS, n=336) to investigate relationships between maternal psychopathology, TVC of socioeconomic status, and offspring depressive symptoms. Results: We found that, on average, offspring depressive symptoms score increased by 3.9% (95% CI: 1.0%-6.9%, p = 0.039) for each unit of maternal psychopathology (SRQ) at 48 months whilst adjusting for time-varying socioeconomic status (at 18, 30, 42 and 54 months). Our simulations identified several realistic scenarios where selections ignoring TVC - with TVC mediated exposure effects present - were prone to be incorrect, including our DCHS example. Conclusion: DME-SLCMA is a robust new approach for life course modelling in the presence of time-varying covariates. We recommend adjusting for TVC whenever possible, and, when not possible, our simulation study identified that scenarios where mediated effects are comparable, or greater, in magnitude to direct effects are most prone to confounding.

14
Multi-ancestry genome-wide association study and meta-analysis of stimulant use disorder reveals biology and relationships to other psychiatric disorders

Beck, S. E.; Deak, J. D.; Levey, D. F.; Ge, T.; Jeffries, P. W.; Lai, D.; Mallard, T. T.; Degenhardt, L.; Lind, P. A.; Tollerup Nielsen, T.; Tubbs, J. D.; Wetherill, L.; Johnson, E. C.; Hatoum, A. S.; The SUD Working Group of the Psychiatric Genomics Consortium, ; COGA Collaborators, ; Yale-Penn Collaboration, ; The VA Million Veteran Program, ; Borglum, A.; Demontis, D.; Medland, S. E.; Martin, N. G.; Nelson, E. C.; Smoller, J. W.; Kranzler, H. R.; Gaziano, J. M.; Stein, M. B.; Agrawal, A.; Edenberg, H. J.; Gelernter, J.

2026-06-10 genetic and genomic medicine 10.64898/2026.06.05.26354997 medRxiv
Top 1%
1.2%
Show abstract

Stimulant use disorder (StimUD) is a significant public health problem, but genetic studies have been limited by small sample sizes. We conducted genome-wide association studies (GWAS) of StimUD in the Million Veteran Program (MVP) and All of Us (AOU), followed by meta-analysis with FinnGen and 10 additional datasets, for a total of 709,369 individuals (Ncases=33,977, Ncontrols=675,392) in four broad ancestry groups: European (EUR) (Ncases=22,564, Ncontrols=624,672), African (AFR) (Ncases=7,574, Ncontrols=34,189), Admixed American (AMR) (Ncases=3,657, Ncontrols=15,698), and East Asian (EAS) (Ncases=182, Ncontrols=833). Population-specific SNP heritability was 6.1% in EUR and 2.4% in AFR. We discovered a total of 19 genome-wide-significant loci, six in EUR, including DRD2*rs5794864, P=7.32E-10, one in AFR, five in a multi-ancestry meta-analysis, including CHRNA5*rs55781567, P=3.27E-9, two in a male-only meta-analysis, including FTO*rs8057044, P=9.50E10-9, and five in a meta-analysis of sex-stratified results. In a hold-out AOU subsample (NEUR=18,841, NAFR=12,263, NAMR=9,739), ancestry-specific polygenic risk scores were significantly associated with StimUD in EUR (OR=3.28, 95% confidence interval (CI)=2.89-3.71) and AMR (OR=2.01, 95% CI=1.71-2.37). Transcriptome-wide association studies, fine-mapping, and colocalization analyses prioritized additional genes (e.g., GPX1, BSN). Genetic correlation, Mendelian randomization, and causal mixture analyses revealed relationships with other substance use and use disorder phenotypes, including cannabis use disorder (rg=0.94, P=5.43E-237) and opioid use disorder (rg=1.01, P=4.40E-107), and other psychiatric traits, including anxiety, depression, neuroticism, and attention-deficit/hyperactivity disorder. This is the first well-powered GWAS of StimUD, and it offers significant insights into disease biology.

15
Metatranscriptomics-Derived Disease Risk Scores as a Preventive, Diagnostic, and Treatment Support Tool

Hu, L.; Bass, M.; Patridge, E.; Molusky, M.; Antoine, G.; Vuyisich, M.; Banavar, G.

2026-06-06 genetic and genomic medicine 10.64898/2026.05.29.26354333 medRxiv
Top 1%
1.1%
Show abstract

Background: Chronic diseases and symptom syndromes often develop after prolonged biological changes that may precede formal diagnosis. RNA-based metatranscriptomics captures active microbial and human gene expression and may provide a functional layer for disease risk evaluation. To address this translational gap, we developed and validated a Disease Risk Score (DRS) framework that integrates metatranscriptome-derived pathway activity scores from stool, saliva, and blood samples, and evaluated its potential clinical utility as an adjunct risk-evaluation tool. Methods: DRS uses disease-specific sets of pathway activity scores derived from stool and saliva microbial functions, stool and saliva microbial taxa, and blood human gene expression. For each disease, 'not optimal' pathway scores are aggregated into a normalized cumulative odds ratio, or cOR, using score-level odds ratios, statistical significance, and literature-supported biological relevance derived from a Development Cohort of 22,369 individuals. A cOR [&ge;] 5 is defined as high risk. Performance is evaluated in an independent Validation Cohort of 15,908 individuals using self-reported diseases as the reference. Disease support requires both significant cOR separation between self-reported and not-reported (Cohen's d [&ge;] 0.2) and risk ratio enrichment of self-reported disease among individuals classified as high risk (95% CI of Risk Ratio > 1). Results: Of 20 initially evaluated diseases, 15 meet the prespecified validation criteria on the independent validation cohort: ADHD, anxiety, chronic fatigue syndrome, depression, GERD, hypertension, inflammatory bowel disease, IBS-C, IBS-D, insomnia, MASLD, obesity, obstructive sleep apnea, Sjogren's syndrome, and type 2 diabetes. Five selected clinical scenarios illustrate how DRS can support clinician-mediated decision making, including IBS subtype reclassification, improved diagnostic acceptance in IBS-D, personalized lifestyle counseling in MASLD and early type 2 diabetes, and diagnostic uncertainty in atypical GERD. Conclusions: DRS is a metatranscriptomics-based risk-stratification framework that aggregates active microbial and human pathway signals into interpretable disease-specific risk estimates across a wide range of disease conditions. Validation against self-reported disease labels in an independent cohort shows significant risk enrichment for each of 15 diseases. DRS is intended as an adjunct to clinical evaluation: a decision support tool in situations where routine care encounters uncertainty, delay, or low patient engagement. Future prospective studies using clinically adjudicated endpoints are needed to assess calibration and clinical outcomes.

16
Identifying Clinical Diagnostic Trajectories Associated With Suicide Death Using Temporal Sequence Mining of Linked Claims and Mortality Data

Belouali, A.; Kitchen, C.; Haroz, E.; Lehmann, H.; Nestadt, P. S.; Wilcox, H. C.; Kharrazi, H.

2026-06-10 health informatics 10.64898/2026.06.08.26355231 medRxiv
Top 1%
1.0%
Show abstract

Background: Most approaches to suicide risk assessment consider clinical conditions as independent risk factors, potentially overlooking prognostic information in the order in which conditions accumulate. We applied temporal sequence mining to linked claims and mortality data to identify ordered clinical diagnostic trajectories associated with suicide death. Results: The cohort included 3 647 059 insured Maryland residents aged 10 years or older with available claims records in the Maryland Suicide Data Warehouse from January 1, 2016, to December 31, 2020, among whom 768 suicide deaths were ascertained through medical examiner linkage. Sequential pattern mining of ICD-10-CM diagnoses grouped into Clinical Classifications Software Refined categories identified 89 221 candidate sequences, of which 1 816 remained significantly associated with suicide death in time-varying Cox models. Adjusted hazard ratios (AHRs) ranged from 2.4 to 134.1. Two-thirds of significant trajectories ended in physical conditions, and approximately half crossed from psychiatric to physical endpoints. Among suicide decedents, 62% were exposed to at least 1 significant sequence (median, 16 per case); median sequence duration was 18.7 months, and median time from completion to death was 13.1 months. In landmark analyses, among patients with depression who later developed suicidal ideation (n = 26 356), the path through anxiety, then anemia, was associated with higher risk (AHR, 4.6; 95% CI, 2.2-9.5), whereas the anxiety-only path was not (AHR, 1.3; 95% CI, 0.8-2.1). Among patients with anxiety who later developed hypertension (n = 149 215), the path through history of self-harm was associated with higher risk (AHR, 32.0; 95% CI, 16.6-61.6). Associations were generally consistent across sex and age. Conclusions: Temporal ordering of clinical conditions may carry prognostic information for suicide death. Clinical trajectories incorporating physical illness within psychiatric sequences identified higher-risk groups. These findings suggest that opportunities for risk detection may extend beyond psychiatric settings and that suicide risk signals may be fragmented across care settings and not apparent within isolated encounters.

17
Inflammation and late life depressive symptoms

Forbes, M.; Lotfaliany, M.; Miteku, B. M.; Yu, C.; Lacaze, P.; Isvoranu, A.-M.; Kang, M.; Nguyen, T.; Woods, R.; McNeil, J.; Neumann, J.; Mohebbi, M.; Berk, M.

2026-06-10 psychiatry and clinical psychology 10.64898/2026.06.05.26354416 medRxiv
Top 1%
0.9%
Show abstract

Background Low-level systemic inflammation has been associated with late-life depressive symptoms. Whether individuals with higher inflammation derive preventive benefit from low-dose aspirin therapy is unknown. Methods We performed a post-hoc analysis of the ASPiring in Reducing Events in the Elderly (ASPREE) randomised, double-blind, placebo-controlled trial. Baseline C-reactive protein (hsCRP) was measured in plasma and depressive symptoms were assessed annually using the Center for Epidemiologic Studies Depression 10 Scale with elevated symptoms defined as CES-D-10 >= 8. Participants with elevated depressive symptoms at baseline were excluded. We fitted population-averaged logistic generalised estimating equation models adjusted for baseline sociodemographic and lifestyle covariates, including an hsCRP x treatment interaction to test effect modification by aspirin. Results Higher baseline hsCRP was associated with increased odds of elevated depressive symptoms during follow-up (OR 1.07 per SD increase in hsCRP, 95% CI 1.03-1.11). Low-dose aspirin allocation did not modify the hsCRP-depressive symptoms association (interaction OR 1.02, 95% CI 0.94-1.10). Findings were similar after additional adjustment for comorbidity and other covariates. Conclusions In community-dwelling older adults during the ASPREE randomised trial period, higher baseline hsCRP was modestly associated with elevated depressive symptoms. There was no evidence that low-dose aspirin was associated with reduced risk of depressive symptoms among participants with higher baseline inflammation.

18
Multivariate Machine Learning Analysis of M-ECG-derived Heart Rate Variability in TBI Veterans, With and Without Comorbid PTSD

Izadysadr, A.; Bagherzadeh, H. S.; Rowland, J.; Martindale, S. L.; Stapleton-Kotloski, J. R.; Godwin, D.

2026-06-08 psychiatry and clinical psychology 10.64898/2026.06.05.26354915 medRxiv
Top 1%
0.8%
Show abstract

Traumatic brain injury (TBI) and posttraumatic stress disorder (PTSD) frequently co-occur in Veterans, producing overlapping symptoms and shared autonomic dysregulation. Heart rate variability (HRV) offers a noninvasive measure of autonomic function. Univariate HRV analyses often fail to capture complex, multivariate patterns associated with comorbidity. This study applied machine learning to HRV features extracted from MEG-derived electrocardiogram (M-ECG) signals to differentiate Veterans with TBI alone (TBI-alone; n = 42) from those with comorbid PTSD (TBI+PTSD; n = 40). Time-domain, frequency-domain, geometric, and nonlinear HRV metrics were analyzed using nested cross-validated Random Forest and XGBoost classifiers, with Boruta-based feature selection and SHapley Additive exPlanations for model interpretability. Both classifiers achieved above-chance discrimination (Random Forest AUC = 0.663; XGBoost AUC = 0.635). Multivariate models identified distributed autonomic signatures in TBI+PTSD, including altered sympathovagal balance, increased low-frequency proportion, and greater heart rate complexity. In contrast, univariate HRV differences were subtle and did not survive correction for multiple comparisons. These findings demonstrate how using multivariate machine learning HRV analysis could help with detecting comorbidity-specific autonomic patterns, suggesting that HRV-derived signatures may serve as exploratory biomarkers for risk assessment and targeted interventions in Veterans with TBI and PTSD.

19
Global and local genetic overlap among ME/CFS, irritable bowel syndrome and psychiatric traits: a hypothesis-generating analysis

Lee, J.

2026-06-10 psychiatry and clinical psychology 10.64898/2026.06.08.26355171 medRxiv
Top 2%
0.7%
Show abstract

Background. Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) and irritable bowel syndrome (IBS) frequently co-occur following infection, yet shared genetic architecture at the locus level has not been systematically characterised. Aims. To estimate global and local genetic correlations between ME/CFS (including infection-onset subgroup), IBS, major depressive disorder (MDD) and loneliness/isolation, and characterise ME/CFS cell-type heritability enrichment. Method. GWAS summary statistics: DecodeME (15,579 ME/CFS; 9,738 infection-onset), FinnGen R9 (9,296 IBS), PGC MDD Wave 2 (45,396) and UK Biobank loneliness (N=455,364). LDSC for global correlations; LAVA for local correlations across 2,495 loci; MAGMA for cell-type enrichment (Descartes Human atlas); coloc.abf for colocalisation. Results. All pairwise global correlations were significant after Bonferroni correction, including ME/CFS-all-MDD (rg=0.598, 95% CI 0.46-0.74) and ME/CFS-all-IBS (rg=0.573, 0.39-0.75). Of 4,232 local tests, 16 reached FDR<0.05; two lonelinessxMDD loci were Bonferroni-significant. ME/CFS-MDD showed three FDR-significant local correlations, but all were boundary-estimated and non-Bonferroni-significant. A borderline infection-onset ME/CFS-IBS signal occurred at chr12q24.22 ({rho}=1.000, FDR=0.046), but colocalisation did not support a shared causal variant (PP.H4=0.007). ME/CFS heritability was enriched in inhibitory neurons (P=1.210x-7) and enteric nervous system neurons (FDR=0.004), with no FDR-significant peripheral immune cell-type enrichment in the atlas used. Conclusions. High global ME/CFS-MDD correlation was accompanied by limited, boundary-estimated, non-Bonferroni-robust local sharing; the data do not support reducing ME/CFS to depression at the genetic-architecture level. Neural enrichment, including enteric nervous system neurons, supports involvement of neural components in ME/CFS susceptibility without excluding immune mechanisms. A borderline ME/CFS-IBS signal at a NOS1-containing region generated hypotheses requiring replication.

20
The polygenic risk score and inter-familial heterogeneity in multigenerational families affected by schizophrenia and bipolar disorder

Ricard, J.; Dubeau, A.; Moreau, C.; Boisvert, M.-C.; Maziade, M.; Bureau, A.; Girard, S. L.

2026-06-08 psychiatry and clinical psychology 10.64898/2026.06.08.26354912 medRxiv
Top 2%
0.7%
Show abstract

In the past two decades, the focus on genome-wide association studies in large samples of unrelated patients has overshadowed family genetic studies. Therefore, little is still known about the levels and effects of the transmission of polygenic risk scores (PRS) among familial cases of schizophrenia (SZ) or bipolar disorder (BD) and their unaffected relatives. Prior research has shown that PRS are elevated in both patients and young individuals at familial risk for BD and SZ. We sought to study the transmission of PRS in affected multigenerational families and non-affected adult relatives (NAARs) with or without other non-mood nonpsychotic DSM-IV diagnoses and unrelated non-affected individuals from the same population. We genotyped 1,117 participants divided in 48 families from the Eastern Quebec Schizophrenia and Bipolar Disorder Kindreds. PRSs for both SZ and BD were computed using Multivariate Lassosum. For both SZ PRS and BD PRS, SZ and BD cases present higher PRS compared to controls, replicating previous findings. Regardless of a diagnosis of other non-psychotic and non-mood conditions, NAARs presented higher PRS than the unrelated cohort. Crucially, a subset of families presented consistently low PRS transmission profiles across generations, falling below expectations from our polygenic inheritance model. When the effect of individual PRs is accounted for, we observed sex-specific associations between familial PRS and patients' symptom dimensions. Our results clearly demonstrate that polygenic inheritance alone does not adequately explain disease transmission in families. Such an approach may also clarify why some families exhibit dense clustering of cases despite minimal polygenic burden.